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The non-topological FL model is studied at finite temperature and density. The soliton solutions 
of the FL model in deconfinement phase transition are solved and thoroughly discussed for different 
boundary conditions. We indicate that the solitons before and after the deconfinement have different 
physical meanings: the soliton before deconfinement represents hadron, while the soliton after the 
deconfinement represents the bound state of quarks which leads to a sQGP phase. The corresponding 
phase diagram is given. 
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I. INTRODUCTION 

An important discovery at Relativistic Heavy Ion Collider (RHIC) in recent years is the realization of strongly 
interacting Quark Gluon Plasma (sQGP) For the collective effects from RHIC experiments, known as radial 

and elliptic flow, the QGP could be well described by ideal hydrodynamics. It implies that the QGP at RHIC is the 
most perfect fluid j^-Q ■ For this reason, Shuryak has pointed out there should be lots of bound states [l| especially for 
light and heavy qq bound states at Tc < T < ATc, where Tc is the transition temperature. For the physical quark mass 
the QCD transition is a non-singular crossover at finite temperature. In recent lattice results, at high temperatures 
and small densities, it is a crossover from hadronic phase to QGP phase 0, Q. From recent theoretical studies, we 
have realized the rich phases of QCD theory, such as color superconducting phase 0, pion condensation color 
glass condensate and quarkyonic phase [ll|- Now the RHIC results tell us that after deconfinement, the quarks are 
not free immediately, as we usually thought for a weakly coupled QGP (wQGP), but still in strong coupled state 
which leads to a new state of nuclear matter as the sQGP. How to understand the formation mechanism of sQGP is 
an important problem fl^ . As we know, there are confine potential and color coulomb potential between quarks in 
vacuum, know as V a/r + kr. With temperature T increased to some critical temperature, the confine potential 
disappears, while the coulomb potential remains, which could still be very strong 131. Thus in literatures, the effective 
coulomb potential has been often used to describe the strong interaction of sQGP 14]. But from the theoretical point 
of view, it is lack of a more fundamental theoretical model to understand these bound states of quarks. In this paper 
we wish to use non-topological FL model to study this problem. 

FL model has been widely discussed in past decays |15l-ll7| . It has been very successful in describing phenomenolog- 
ically the static properties of hadrons and their behaviors at low energy. The model consists of quark fields interacting 
with a phenomenological scalar field a. The a field is introduced to describe the complicated nonperturbative features 
of QCD vacuum. It naturally gives a color confinement mechanism in QCD theor y. T he model has been also extended 
to finite temperatures and densities to study deconfinement phase transition |18l421| . However it seems that the deep 
meaning of the soliton solutions in deconfinement phase transition has not been revealed in the past studies. The main 
purpose of this paper is to study in detail the properties of the solitons in FL model before and after deconfinement, 
and provide a natural explanation and description of sQGP through this effective theoretical model. 

The organization of this paper is as follows: in section 2 we give a brief introduction of FL model. The field 
equations and the effective potential are derived. In section 3, the soliton equation of FL model is solved for different 
boundary conditions and the physical meanings of these solitons in deconfinement are thoroughly discussed. In section 
4, a phase diagram of deconfinement phase transition is given. The last section is the summary. 



II. FL MODEL AND THE EQUATIONS 



We start from the Lagrangian of the FL model. 
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'0 represents the quark field, and a denotes the phenomenological scalar field, a, b, c, g and B are the constants which 
are generally fitted in with producing the properties of hadrons appropriately. 

For the thermal equilibrium system, the a field will be replaced by ct + cr', where a is the Gibbs thermal average of 
the a field, and a' is the fluctuation. At mean field approximation, one can obtain the field equations as follows, 

(Z7^9'^ - ga)i^ = 0, (3) 
dV.ff 



da 
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where V^f / is one loop effective potential at finite temperature and density. Notice that at mean field approximation the 
quantum corrections have been neglected. The V^ff could be derived through finite temperature field theory [20l [2l|. 
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in which [3 is the inverse temperature, is the chemical potential, and 7 is a degenerate factor, 7 — 2{spin) x 
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2{ flavor) x 3{color). Besides Ec — \/ + m'^, Eq — ^ jp + rUq, where ruq = ga and = 

effective masses of the quark and a field respectively. 

From equation Q , one could see that the properties of the soliton field a depend completely on the effective potential 
Vef / . As known, by the effective potential at finite temperature and density, one can study the deconfinement phase 
transition in FL model. Thus the properties of solitons in deconfinement, especially the relations between solitons 
and deconfinement, could be well studied through solving the equation (|4]) at finite temperature and density. 

In our calculation, the parameters are chosen to be a = 17.7 fm~^,b = — 1457.4/m~^, c = 20000,(7 = 12.16. The 
effective mass of a field is fixed at mg. = 550MeV [l^l. At zero temperature and density, the 14// is just the U{(f) 
and could be plotted in FigH] There are two minima: one corresponds to the perturbative vacuum at = 0, another 
corresponds to the physical vacuum at a = a^. 
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FIG. 1: The effective potential at zero temperature and density . 



III. THE SOLITON SOLUTIONS AND SQGP IN DECONFINEMENT 

For the static and spherically symmetric soliton field, the equation ([4]) becomes 

d^a 2 da dVeff 

1 = — (6) 

dr-^ r dr da 

We will solve this equation numerically. A standard numerical package COLSYS will be used for the numerical 
calculation [2^ . 

For solving the soliton equation ([B]), one should first determine the boundary condition. There are two different 
cases of the boundary conditions for solving the equation [l3, H^. In the following the two cases will be studied 
separately. 
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A. Boundary condition I 



In this case, we may consider the mechanical analog in which there is a point particle at "positon" a and a time "r" , 
moving in a "potential" —Vefj. The discussions about this mechanical analog could be found in Ref. [23| . Equation 
is then the movement equation of the particle. Because the energy conservation will be applied in this case, the 
second term at l.h.s of equation ([6|) will be omitted. Thus the equation becomes 
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da 



(7) 



In order to determine the boundary condition, suppose at "time" r — —oo, the "position" of the particle is at = 0. 
From Figlll if the particle is pushed very gently to start moving, the particle will move to another point just at 
a = aa where V'e//(0) = Veff{(Ta), noticing that the particle is moving in the "potential" —Veff. Because of energy 
conservation, the particle will then move back and return to the point a — 0. Thus one could obtain the following 
boundary condition. 



cr(r = — oo) = 0, (T(r = +oo) = 0. 



(8) 



From FiglU the physical vacuum a = is stable, and B — Ve//(0) — ^effio'v) is the bag constant. The quarks are 
confined in a soliton bag. At zero temperature and density, the soliton solution represents a hadron. 
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FIG. 2: The effective potential for different temperatures: Ti = 50MeV, T2 = lOOMeV, Td = 121MeV, T3 = ISOMeV, 
T4 = UOMeV and Tc2 = USMeV. 

At finite temperature, the Veff could be plotted in Figl2] There are two critical temperatures Td and Tc2- At Td, 
the two minima degenerate and the bag constant B = 0. At Tc2, the minimum tr = (t„ just vanishes. For different 
configurations of Ve//, the equation ([7]) could be numerical solved by the numerical package COLSYS. However the 
Veff is highly non-linear in a as shown in equation ([S]). It could not be directly used in COLSYS evaluation. At finite 
temperature, as the Veff could be plotted as shown in FiglJl we can fit the curve in terms of cr^, and a'* with the 
effective parameters a{T), b{T) and c(T). That means 

Veff = ^a{T)a^ + ^biT)a' + ^ciT)a* + B{T). (9) 

In this form, the Veff could be applied in COLSYS for the numerical calculation. In the following we will solve 
equation (O and discuss the soliton solutions for different configurations of Vef / . 

At T < Tci, the physical vacuum at ct = ct^, is stable. The bag constant B and its value decreases with 
temperature increasing. The corresponding quarks are confined. The equation ([7]) could be numerically solved and 
the soliton solutions are plotted in Figl3^. One could see that the solitons at temperature T < Td do not change so 
much when compared to that at zero temperature except the peaks of the solitons become higher when temperature 
increases. The quarks keep confined in a hadronic state until the critical temperature T = Td the two vacuums 
degenerate. At this time, still required by the energy conservation the boundary condition is 



cr(r = —00) = 0, a{r — +00) = a, 



(10) 
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The soliton becomes topological as shown in Figl3^. As the bag constant _B = at this time, the hadrons are 
destructed and the quarks are deconfined. 

If we keep increasing temperature to < T < Tc2, from Figl^J the perturbative vacuum tr = is stable. The 
equation (O will be solved under the boundary condition 



a{r — — oo) 



(11) 



The soliton solutions could be plotted in FiglSb- The soliton solutions still exist but become very different from those 
at T < Tel. The peaks of the solitons become downward. When temperature increasing, they become flat. These 
solitons do not represent hadrons any more but the bound states of quarks. Though the quarks are deconfined, they 
are still in strong coupled states. That is to say the system is in a sQGP phase. Until T > Tc2, the physical vacuum 
a = a.-a vanishes and only the perturbative vacuum = exists, has the soliton solution disappeared. At this time 
the quarks become free. The system goes into a quasi-free gas phase of QGP. 
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FIG. 4: The effective potential for different chemical potentials: /ii = lOOMeV, \i2 = 200MeV, ^ici = 255MeV, [iz = SOOMel/, 
= SSOMeV and ^,,2 = 375MeV. 

At finite temperature and density, we could also make the similar discussions. In FigUl we show the 14// at 
T = h^MeV and different chemical potentials. There are also two critical chemicals /id and [i,c2- At [i^x the two 
vacuums degenerate; at /ic2 the vacuum of cr 7^ just vanishes. The corresponding soliton solutions are plotted in 
FigEK and FiglSb. It is clear that at /i < /Xci, the system is in a hadronic phase; at /id < \x < /ic2, the system is in a 
sQGP phase; at /< > /ic2, the system becomes quasi-free gas of QGP. 
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FIG. 5: The soliton solutions for boundary condition I at different chemical potentials, (a) fii = lOOAdeV, /12 ~ 200MeV and 
Hci = 255MeV. (b) ^3 = 300A//eV and fn = 350MeV. 
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FIG. 6: The soliton solutions for boundary condition II at different temperatures. (a)ro = OMeV, Ti = 50MeV, T2 = WOMeV 
and Tci = UlMeV. (b) T3 = ISOMeV and T4 = UOMeV. 



B. Boundary condition II 



In this case, at T < Td, the boundary condition is taken as 

-0, aU^=0. (12) 

The equation Q could be numerically solved at finite temperature and density for different configurations of 14/ / . 

From FiglH at T < Td , the soliton solutions are plotted in FiglS^. The physical vacuum ct = ct^ is stable and The 
bag constant B ^ 0. The solitons still represent that the quarks are confined and the system is in a hadronic phase. 
Until T — Tci the system are deconfined. The corresponding soliton is also plotted in FiglSK- 

At Tci < T < Tc2, the boundary condition is 

= 0, ijl^^ - a„. (13) 

The soliton solutions are solved and plotted in FiglBlD. These solitons are quite different from those at T < Td- As 
the system is already deconfined, these solitons represent the quarks are in a bound state. The system is in a sQGP 
phase. At r > Tc2, there are no soliton solution any more, the system goes to a quasi- free gas of QGP. 

For different chemical potentials at fixed temperature T — bOMeV, according to the Veff as shown in FiglU the 
corresponding soliton solutions could be solved and plotted in Figl?^ and FiglTb. It is clear that at ^ < fici it is a 
hadronic phase; at ^ci < A* < ^02, it is in a sQGP phase; at /i > it is a quasi-free gas of QGP. 

From above discussions, it is clear that for the two different boundary conditions, the physical meaning of the 
solitons are the same: before deconfinement the solitons represent hadrons while after deconfinment the solitons 
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FIG. 7: The soliton solutions for the boundary condition II at different chemical potentials, (a) = lOQMeV , fi2 = 200MeV 
and Aid = 255MeV. (b) fis = SOOMeV and ^4 = 350MeV. 



represent bound states of quarks. For the two cases, the phase structure of deconfinement remains unchanged and 
unique. 



IV. THE PHASE DIAGRAM 



Now we are in a position to plot the phase diagram of deconfinement phase transition in FL model. When the 
chemical potential is fixed, one could obtain two critical temperatures, fncreasing the chemical potential to another 
fixed value, the other two corresponding critical temperatures could be obtained, and so on. The fi — T phase diagram 
could be plotted as shown in FiglH The full line is the dividing line between hadronic phase and deconfined quark 
phase. The dashed line further divides the deconfined quark phase into the sQGP phase and the wQGP phase. 
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FIG. 8: The phase diagram of deconfinement in the FL sohton model. 



V. SUMMARY 



In this paper the soliton solutions in FL model at finite temperatures and densities have been thoroughly discussed 
for two kinds of different boundary conditions. The physical interpretations of the solitons in deconfinement phase 
transition have been given. As in FL model the system is deconfined at the time that the bag constant becomes zero, 
we indicate that the solitons before deconfinement are hardrons, while the solitons after deconfinement represent the 
strong coupled bound states of quarks, which leads to the sQGP phase. The QCD phase structure still needs lots of 
thorough investigations, especially for the finite density areas, where the results are often model dependent. Here we 
only present a possible phase diagram based on FL model. The phase structure of deconfinement phase transition 
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here is qualitatively consistent with that of Shuryak The whole phase diagram is divided into three phases: the 
hadronic phase, the sQGP phase and the wQGP phase. 
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